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Abstract 

A comprehensive theory of the Weyl-Wigner formalism for the canonical pair angle-angular 
momentum is presented. Special attention is paid to the problems linked to rotational periodicity 
and angular-momentum discreteness. 



1. Introduction 

Phase-space methods were introduced in the very early times of quantum mechanics to avoid 
some of the troubles arising in the abstract Hilbert-space formulation. The pioneering works 
of Weyl [1|, Wigner [2|, and Moyal [3| paved thus the way to formally representing quantum 
mechanics as a statistical theory on phase space Bll5ll6ll7l[8ll9l [T0l[TT1l . 

In this approach one looks for a mapping relating operators (in Hilbert space) to functions (in 
phase space). In particular, quasiprobability distributions are the functions connected with the 
density operator ||T2] [13] [14] [T5] [T%1 [T7J [T8]l . Usually, the relevant theoretical tools are illustrated 
with continuous variables (such as Cartesian position and momentum), for which much physical 
knowledge has been inferred. Likewise, the increasing role of qudits in modern quantum infor- 
mation has fuelled a lot of interest in consistently extending these achievements to the case of a 
discrete phase space fl^l^l2Ul^l^l2^l25lE6ll27ll28ll29l. 

However, other physical problems call for different topologies of the phase space [ 30, 3T1I321 
[33l . Specifically, in this paper we focus on the discrete cylinder Si x Z (Si denotes here the unit 
circle and Z the integers), which is associated to the canonical pair angle-angular momentum. 
A number of systems, such as molecular rotations, electron wave packets, Hall fluids, and light 
fields, to cite only a few examples, can be described in terms of this geometry 11341 l35l . In 
quantum optics, it is of primary importance to deal with the orbital angular momentum of twisted 
photons ||36l[37l , which have been proposed for numerous applications ||38ll39ll . 

The proper definition of angles in quantum mechanics has a long history and requires more 
care than perhaps might be expected B01I4TI I42 43 1. Since the conjugate angular momentum has 
an unbounded spectrum (that includes positive and negative integers), it is, in principle, possible 
to introduce a bona fide angle operator. Periodicity, however, brings out subtleties that have 
triggered long and heated discussions. We think that a phase-space treatment of this variable 
may shed light on the origin of all those problems. 

A pioneer attempt in that direction was made by Mukunda Il44l l45ll . who for the first time 
constructed a Wigner function on the discrete cylinder. This work was subsequently reelaborated 
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and developed in a variety of directions by other authors ||46l l47l [48 , 49 ,50 5 1 ] . Coherent states 
for this pair have also attracted much attention 021 l53l ,54 55 , 56 . 57 1 . These accomplishments, 
as well as some other related questions, such as the associated minimum uncertainty states, are 
reviewed in great detail by Kastrup [|58ll . 

In spite of this considerable progress, there is still the need, in our view, for a comprehensive 
geometrical approach to this problem, which allows us to integrate in a structured way all the 
previous knowledge, while providing the right tools to attack recent and attractive problems, 
such as, e.g., the tomography of these systems |59l . The purpose of this paper is precisely to 
provide such a description. 

2. Phase space for quantum continuous variables 

In this section we recall the structures needed to set up a phase-space description of Cartesian 
quantum mechanics. This will facilitate comparison with the cylinder later on. For simplicity, 
we choose one degree of freedom, so the associated phase space is the plane M 2 . 

The relevant observables are the Hermitian coordinate and momentum operators q and p, 
with canonical commutation relation (with h = 1 throughout) 

[q,p]=ii, (1) 

so that they are the generators of the Heisenberg-Weyl algebra [60|. Ubiquitous and profound, 
this algebra has become the hallmark of noncommutativity in quantum theory. 

To avoid technical problems with the unboundedness of q and p, it is convenient to work with 
their unitary counterparts 

U (q) = e\p(-iq p) , V(p) = exp(i> q) , (2) 
whose action in the bases of eigenvectors of position and momentum is 

U(q')\q) = \q + q), V(p')\p) = \p + p) , 0) 

so they represent displacements along the corresponding coordinate axes. The commutation 
relations are then expressed in the Weyl form ||6T1 

V(p)U(q) = e'w U(q)V(p) . (4) 

Their infinitesimal version immediately gives ([T]), but Q is more useful in many instances. 
In terms of U and V a general displacement operator can be introduced as 

D(q, p) = U{p)V{q)e- >qpl1 = exp[i(pq - qp)] , (5) 

with the parameters (q, p) labeling phase-space points. These operators form a complete trace- 
orthonormal set (in the continuum sense) in the space of operators acting on H (the Hilbert space 
of square integrable functions on R): 

TAb(q,p)b\q',p')] = 2n6{q - q')6(p - p') . (6) 

Note that fy^{q,p) = D(-q, -p), while D(Q, 0) = 1. In addition, they obey the simple composi- 
tion law 

D(q, p) b(q',p') = exp[i(qp' - pq')/2] D(q + q' , p + p') . (7) 
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We can also work with the Fourier transform of D(q, p) 

,P) = f exp[-i(pq' - qp')]D(q',p')dq'dp' , (8) 

(2?r) 2 J R 2 

which is called a Stratonovich-Weyl quantizer [4|. One can check that the operators w(q,p) are 
also a complete trace-orthonormal set that transforms properly under displacements 



where 



and 



w(q,p) = D(q,p)w(0,0)DHq,p), (9) 

w(0, 0) = I D(q, p)dqdp = 2P, (1 0) 

Jr 2 

P= f \q)(-q\dq= f \p)(-p\dp (11) 
Jr Jr 

is the parity operator. 

Let A be an arbitrary (Hilbert-Schmidt) operator acting on ti. Using the Stratonovich-Weyl 
quantizer we can associate to A a tempered distribution a(q, p) on M. 2 representing the action of 
the corresponding dynamical variable in phase space. In fact, this is known as the Wigner-Weyl 
map and reads as 

a(q,p) = TT[Aw(q,p)]. (12) 

The function a(q, p) is the symbol of the operator A. Conversely, we can reconstruct the operator 
from its symbol through 

A = — —r a(q,p)w(q,p)dqdp . (13) 
(2tt) 2 J R 2 

The symbols form an associative algebra endowed with a noncommutative (star) product 
inherited by the operator product. Let the functions a(q, p) and b(q, p) correspond to the operators 
A and B, respectively. If we denote by {a*b)(q, p) the function corresponding to AB, some simple 
manipulations lead us to 1 16 1 

(a*b)(q,p)- r- f f a(q+q' , p+p') exp[2i(q' p" -q" p')]b(q+q" , p+p") dq'dp'dq"dp" . 

(14) 

Expanding the functions a and b in a formal Taylor series at the point (q, p) e K 2 we obtain 

(a*b)(q,p) = a(q,p) exp(-^) b(q, p) , (15) 

<— > 

where f is the Poisson operator 

dq dp dp dq ' 
in terms of which we can define the Moyal bracket J3] 

\a,b}M — -(a ★ b — b ★ a) , (17) 
i 



which is the phase-space counterpart of the commutator in the Hilbert-space formulation. Con- 
sequently, the time evolution of a(q, p) can be expressed as 

da 

- = (M1«, (18) 
at 

where h is the symbol of the Hamiltonian. 

In this context, the Wigner function is nothing but the symbol of the density matrix g. There- 
fore, 

Wg(q,p) = Tr[gw(q,p)] , 

(19) 



In) 1 X 



(2. 

For a pure state \ X Y), it can be represented as 



e=7^y I w(q,p)Wg(q,p)dqdp . 



W e (q,p) = ^- f exp(ipq')V(q-q'/2)y(q + q'/2)dq' , (20) 



which is, perhaps, the most convenient form for actual calculations. According to ( 18 i, the 
Liouville-von Neumann equation in phase space is 

dWg 

-^-={h,W e ] M . (21) 

When applying this equation in practical cases, one needs the following mapping for the action 
of the basic variables q and p on g: 

qQ >-> [ q ~^~^ W e^P^ PC [p + ~2~^) W e ( ' q ' l ^ ' 



Qq >-> (q+ \^~\ w M>p)> eP ^ [ p ~ ^tq) Ws< ~ q ' l ^ ' 



(22) 



The Wigner function defined in ( [T9| fulfills all the basic properties required for any good 
probabilistic description. First, due to the Hermiticity of w(q,p), it is real for Hermitian op- 
erators. Second, on integrating W(q,p) over the lines q^ = qcosfi + psinfi, the probability 
distributions of the rotated quadratures q$ are reproduced 

I W e (q, p) 5(q - q$) dqdp = (q^glq^) ■ (23) 
Jr 2 

In particular, the probability distributions for the canonical variables can be obtained as the 
marginals 

f W e {q,p)dp = {q\Q\q), f W e (q, p) dq = (p\g\p) . (24) 
Jr Jr 

Third, W e (q, p) is translationally covariant, which means that for the displaced state g' = D(q',p') 
g D^iq' , p'), one has 

Wty(q,p) = W s (q -4,p-p'), (25) 
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so that it follows displacements rigidly without changing its form, reflecting the fact that physics 
should not depend on a certain choice of the origin. The same holds true for any linear canonical 
transformation. 

Finally, the overlap of two density operators is proportional to the integral of the associated 
Wigner functions: 

Ti(gg)cc I Wg(q,p)Wg>(q,p)dqdp. (26) 

Jr 2 

This property (known as traciality) offers practical advantages, since it allows one to predict the 
statistics of any outcome, once the Wigner function of the measured state is known. 

However, the Wigner function can take on negative values, a property which distinguishes 
it from a true probability distribution. Indeed, this negativity is associated with the existence of 
quantum interference, which itself may be identified as a signal of nonclassical behavior (62). 
The characterization of quantum states that are classical, in the sense of giving rise to nonnegative 
Wigner functions, is a topic of undoubted interest. Among pure states, it was proven by Hud- 
son [63 1 that the only states that have nonnegative Wigner functions are Gaussian states l64ll65l . 

Coherent states are closely linked with the notion of Gaussian states. The displacements con- 
stitute a basic ingredient for their definition: indeed, if we choose a fixed normalized reference 
state IT'o), we have ll66l 

\q,p) = D(q,p)Wo), (27) 

so they are parametrized by phase-space points. These states have a number of remarkable 
properties inherited from those of D(q, p). In particular, D(q, p) transforms any coherent state in 
another coherent state: 

D(q, p) \q, p) = exp[i(q'p - p'q)/2] \q + q',p + p) . (28) 

The standard choice for the fiducial vector I^q) is the vacuum |0). This has quite a number of 
relevant properties, but the one we want to stress for what follows is that |0) is an eigenstate of 
the Fourier transform (as they are all the Fock states) [67 , 68 1. In fact, over this apparently trivial 
property rests a huge amount of physical knowledge. So, |*Po) is taken as the Gaussian 

^o(4)= 4/4 ex P(V/2), (29) 

in appropriate units. In addition, this wave function, as any coherent state, represents a minimum 
uncertainty state, namely 

(Aq) 2 (Ap) 2 = 1, (30) 

where {Aq) 2 and (Ap) 2 are the corresponding variances. 

As a last illustration of the role played by the displacement operators, we note that their 
completeness allows us to write for any observable A 

A = ±- f Tv[ADHq,p)]D(q,p)dqdp. (31) 
27T Jr2 

Evaluating the trace in terms of the set of eigenvectors of the rotated quadratures q^ we obtain l69l 

Wg(q,p)=7T- I d<p I p(qd>,<f>) K(q cos (p +psin</>- q^) dq^dip , (32) 
27T Jo J_ OT 



where K(x) is the kernel 

K(x)= ^limRe 1 (33) 

2 i;|0 (x + IT]) 1 

The tomograms p(q,p, </>) = (q<p\§\q#) are measured in a homodyne detector. This shows that these 
rotated quadratures provide a complete quorum to reconstruct the density operator and hence the 
Wigner function. 



3. Phase space for angle-angular momentum 

3.1. Geometrical properties of the discrete cylinder 

In this Section we trace the changes required when we replace the Cartesian position by 
a periodic angular position <p e S\. For definiteness, we take the window (—n,n) in further 
considerations. The canonical conjugate variable, denoted now as L, is the component of the 
angular momentum along the axis orthogonal to the rotation plane. While classically a point 
particle is necessarily located at a single value of the angle <p, the corresponding quantum wave 
function is an object extended around S\ and can be directly affected by the nontrivial topology. 

One may be tempted to think that angular position should stand in the same relationship to 
angular momentum as ordinary position stands to linear momentum. This would prompt to use 
the commutation relation [0, L] = i 1 and interpret the angle operator as multiplication by <j>, 
while L is the differential operator L = -id$. Nevertheless, the use of this operator may entail 
many pitfalls for the unwary and needs a very subtle analysis ll40l ITTl l72l . 

As stressed in the Introduction, a long discussion about the properties of this angle operator 
turns out to be unnecessary to settle a phase-space description. Indeed, let us denote by \C) 
the basis of angular-momentum eigenstates. In principle, the corresponding eigenvalues can be 
arbitrary. However, when some invariance principles hold [52, 58], we can restrict our attention 
to integer values. In this case, the states 

i0 = _Ljy«*io, 04) 

V2^ 

constitute a basis of states with well-defined angle, so that (<f>\tp') = <52k(<P - 0'X where 
represents the periodic delta function (or Dirac comb) of period 2n Il73l . and they allow for a 
resolution of the identity. Denoting the wave function components in both bases as T^ = {(\V) 
and ¥(0) = (0|Y>, we have 

= -L J] e ^ Ve> = f *e-** dcf, , (35) 

V27T feZ J-n 

which translates the Fourier relationship inherent to any canonical pair. 

To proceed further, we restrict ourselves to the exponentiated version of the canonical pair 

U(4>) = exp(-i</> L) , V(£) = exp(z'€ 0) , (36) 

which, in addition, are experimentally feasible operations lP74ll75ll . We have the action 

0(4>')\4>) = \4> + 4>'), V(t'M = \{ + t'), (37) 
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so that they represent displacements along the coordinate axes of the cylinder Si xZ. Throughout, 
the angle addition and subtraction must be understood modulo 2n. In dealing with this pair, it is 
customary also to take as the fundamental angular variable the complex exponential of the angle 

E = exp(-$) , (38) 

instead of the angle itself. 

The Weyl form for this pair reads as 

WW(<f>) = e u * U(<f>)V(i) , (39) 

and following the ideas of Section|2] the displacement operators can be introduced as 

D{€, cf>) = exp[ia(€, <f>)] V(£)U(4>) , (40) 

where a(l, (p) is a phase. The reasonable condition <p) = £>(—(, -<p) imposes 

exp[ia(£, cp) + ia{-£, -<p)] = exp(-i£0) . (41) 



We cannot rewrite now equation (40 1 as an entangled exponential. 

These displacement operators form a complete trace-orthonormal set 

Tr[L>(A (/>) D\l', cf>')] = 2n5 (e , 5^{<t> - </>') , (42) 

whose resemblance with the relation (|6]l is evident. One can also introduce the Stratonovich-Weyl 
quantizer as 

1 C n 

w(£, <f>) = ——j V exp[-/(^> - ty)] £>(£', </>') d<j>' , (43) 

so that it fulfills the covariance condition 

w(£, <t>) = £>(€, 4>) w(0, 0) £>HA 4>) ■ (44) 
Note carefully that vv(0, 0) cannot be identified now with the parity operator on the cylinder 

£=x/><-^= r i^x-^w. (45) 

as it is the case for continuous variables. This is due to the fact that D(£,<p) P £fi (£,</>) do not 
constitute an operator basis. One gets a basis if one supplements P by EP (or PE), E being the 
exponential of the angle d38l. 

The symbol of an operator A turns out to be 

a{€,<t>) = Tr[Aw({, ((>)], (46) 

while the inversion is given by 



A = 2tt I a(€,<f>)ti/(.€,(f>)d<f>. (47) 
7 



As explained by Plebariski and coworkers [76], this Wigner-Weyl correspondence on the cylinder 
directly distinguishes between the "classical" phase space Si xM from the "quantum" phase space 
Si x Z: only the former admits a formalism based on the use of symbols. 

We can also look for the form of the star product. In close analogy with equation ( |14) , we 
have the integral form 

1 C n 

(a*b)(t,cf>) = — J] I a(e+f, 0+0'/2) exv[i(t"cf>' -£'<(>")] b(M" ,c[>+cf>" /2) d<t>'d<p" . (48) 

Expanding the functions a and b in a formal Taylor series at {(, 0), this result can also be rewritten 
in a "differential" form [compare with equation ( fl5] >] 

(a ★ &)(£, 0) = a(£, 0) exp (-^P* ) 0) . (49) 

where P is the counterpart of the Poisson operator |l6) l on the cylinder 

— > <— 
<— > <— 3 3 — > 

P=5 f — -— 5^. (50) 

O0 O0 

Here d( is some sort of a continuous interpolation between discrete translations, and is defined 
through 

exp(A5t)f{€) = f exp[rX€ + i)0]/(0)^0, (51) 



/(0) being the Fourier transform of the function f(f). 

With this star product, we can introduce a Moyal bracket as in equation ([17} and the time 



evolution of dynamical variables is also given by ( 18 I. 

To conclude, some remarks concerning the phase a((,(f>) in (40 1 are in order. Although 
many of the results in this paper are independent of a{(, 0), a judicious choice of this factor 
may greatly facilitate calculations. The periodicity in requires D((, 0) = D{(, <p + 2n) [apart 
from pT)]. Obviously, there are many admissible functions. For example, quadratic definitions 
like exp[/a'(^, 0)] = exp[-i{(£ + 1)0/2] or exp[/^(^ - 1)0/2] are legitimate, but we discard them 
because they introduce a lot of difficulties when computing practical Wigner functions. 

Another possibility is 

explMA 0)] = exp(-^0/2) exp{/0 [1 - (-l/]/4) , (52) 
leading to the Wigner kernel 

w({, 0) = ^-D({, 0) (P + E P) t)\(, 0) , (53) 
2n 

which, as discussed before, contains the complete set P and EP. In addition, this option brings 
the advantage of very easy calculations, especially for states given in the angular momentum 



representation. However, in this case the term vt>(0, 0) oc P+E P in Eq. ( 53 1 is not parity invariant. 
Since the phase space should not have a "preferred direction", we will also discard this definition. 
We finally fix the phase a as 

exp[ia(€, 0)] = exp(-zT0/2) , (54) 



which may lead to slightly more involved calculations than ( 52 1, but has the clear advantage of 
a parity-invariant kernel. We stress that this definition of a is related to the choice of the 2n 
window for 0, so it is exclusively valid for e [—n, n]. 



3.2. Wigner function on the cylinder 

Since the Wigner function is the symbol associated with the density matrix, we have the 
mapping 

W l >({,cf>) = Tr[ ] Qw({,<f>)], 

(55) 



q = 2tt^ I w(£,(f>)Wg({,</>)d4>. 

fez 

To obtain an operational form of this Wigner function, we note that with the phase definition 
(54 1, the quantizer kernel ( |4"3j ) becomes 



y n )2 X I ( 



WA) = ^ ), f exp[i(^ ' /2 -■£")</>'] ex V [W -t'<l>W"W ~W . (56) 



Because of the contribution with a half-integer exponent, to sum over €' one has to split this sum 
in two different parts according {' is either even or odd. After some manipulations, one gets 

W s ({,cf>) = ±£exp(-2tf»tf~r|£K + 

71 Pel, 

+ ^ X jfpTTn ex P[~'(2^ + l)cf>] H + 1" - t'W + Z" + C + (57) 

2n +1/2 

This expression contains exclusively angular-momentum matrix elements of the density operator 
and it is involved in practical computations due to the second line, which presents a logarith- 
mically slow convergence due to the fact that for odd I, the phase ( 54 1 as a function of <p has a 
discontinuity at <p — +n. Incidentally, this second line would be greatly simplified by employing 
the alternative kernel ( |53| ). 

Sometimes, it is preferable to work in the angle representation, which can be obtained by 
introducing a resolution of the unity in equation (43 i. In this way, we get 

W 6 {€, <t>)=^J exp(/^')<0 - </>'/2\q |0 + </>'/2) d<f>' , (58) 

which coincides with the form introduced in the pioneering work by Mukunda [44|. This Wigner 
function reproduces the proper marginal distributions 

JW« = »), f W^,0)# = <W>, (59) 

and it is explicitly covariant under displacements on the cylinder. This confirms the good prop- 
erties of this approach. 

Note, in closing, that the time evolution of this Wigner function is given by a Liouville-von 
Neumann equation analogous to <\21\ . Now, in practical calculations we shall need the following 
map, which can be easily checked 



(60) 



Again, we have taken as the fundamental angular variable the complex exponential of the angle 
instead of the angle itself. 



4. Examples 

To further appreciate these ideas, we present a few relevant examples. For an angular mo- 
mentum eigenstate and for an angle eigenstate |0o), one has 

W M (€, cf>) = , W m (€, 0) = — 6 27 ,(<f> - 0o) • (61) 

In the first case, the Wigner function is flat in <p and the integral over the whole phase space 
gives the unity, reflecting the normalization of \-£q), while in the second, it is flat in the conjugate 
variable {, and thus, the integral over the whole phase space diverges, which is a consequence of 
the fact that the state \<po) is not normalizable. 

Next, we consider coherent states. The standard approach due to Perelomov |66| does not 
work for the cylinder. However, much as in ( |27| , we can introduce them as 

K0> = />(/,0PPo>. (62) 
The choice of the fiducial vector is not so obvious. We take [561 



i 0s (IW 

¥o(fl = — -r^-A (63) 

where #3 denotes the third Jacobi theta function |77], which incidentally plays the role of the 
Gaussian in circular statistics. In this way, the states \l, <p) inherit the properties of D{£, <p) and 
they turn out to be equivalent to the ones defined in [ 52 1 by an eigenvalue equation or in ll53l via 
a Zak transform. 

The Wigner function for the coherent state |€q,0o) splits as W\e ^y(€, <p) = M,^>) + 



\la,*pa) 

0)- 'Th e " even " turns out to be 



Wg> 0) = l —-e-^ 2 & 3 L - J -) ■ (64) 

This seems a sensible result, since it is a discrete Gaussian in the variable 6, and the equivalent 
for a Gaussian for the continuous angle <p it. However, the "odd" contribution spoils this simple 
picture: 

r_) exp[/(«* -<*o)- 1 /2] / 1 1\ v-i - M exp(-€' 2 - £') 

W f \ At, <p) = 7, , fl 3 {</> ~ <Po + J/2 - > (-1) -r-^ „ , ' ■ 

l4W 2^(0|i) V f +<b-^+l/2 

In figure [TJ the Wigner function for the coherent state \(q = 0, (f>o - 0) is plotted. A pro- 
nounced peak at <p = for t — and slightly smaller ones for £ = ±1 can be observed. The 
associated marginal distributions are also plotted. They are strictly positive, as correspond to 
true probability distributions. For quantitative comparisons, however, sometimes it may be con- 
venient to "cut" this cylindrical plot along a line 0=constant and unwrap it. This is also shown 
in figure [T] 
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Figure 1 : Plot of the Wigner function for a coherent state with = and tpo = 0. The cylinder extends vertically from 
I = -4 to I = +4. The two corresponding marginal distributions are shown. In the right panel we show the unwrapped 
version of the Wigner function. 



A closer look at this figure reveals also a remarkable fact: for values close to <f> — +n and 
( = ±1, the Wigner function takes negative values. Actually, a numeric analysis suggests the 
existence of negativities close to = ±n for any odd value of I. 

As our last example, we address the superposition 

|T> = ^(|A> + ^°K 2 » (65) 
V2 

of two angular-momentum eigenstates with a relative phase e'*>. The Wigner function splits 
again; now the "even" part reads as 

Wjg(*. <f>) = ^{c5f/, + Stjk + 2d (l+ha{ cos[0 o + (h - • (66) 



For the "odd" part, the diagonal contributions vanish, and one has 

Wgto, (f>) = -j cos[0 o + (h - (i)<l>\ =7-*« I+ fe=odd . (67) 



where de l +e 2 =odd indicates that the sum is nonzero only when C\ + 1% is odd. 

In consequence, when \t\ - €%\ is odd, the interference term contains contributions for any I, 



damped as \jt. When \l\ - lj\ is an even number, the contribution (67i vanishes and we have 
three contributions: two flat slices coming from the states and \l->) and an interference term 
located at t= (€ x +{ 2 )/2- 

These features are nicely illustrated in figure[2] The state |*P) is plotted for ti = -3 and I \ = 3 
and and l 2 = -3 and t\ = 4. Changing the relative phase </>q results in a global rotation of the 
cylinder. In can be observed in that the two rings at t = -3 and £ — 4 (as opposed to the rings at 
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I = ±3), are not flat in <f>, but show a weak dependence on the angle due to the odd contributions 
added to the flat Kronecker deltas. 

In view of these results, one can wonder what are the pure state for which the Wigner function 
is nonnegative. The answer was found quite recently [78 1: the Wigner function Wfr>(i, (/)) is 
nonnegative if and only if \*¥) is an angular momentum eigenstate |£ ). The proof is rather 
technical and we skip it. The important point we wish to stress is that while for the continuous 
case the notions of coherent states, Gaussian wave packets, and states with nonnegative Wigner 
functions are completely equivalent, special care must be paid in extending these ideas to other 
physical systems, like angle-angular momentum, since they lose their equivalence. 

To round up this section, we illustrate the benefits of our phase-space approach by dealing 
with the evolution of the "quantum" pendulum 

tf=y + ^(£ + £ 1 ")- (68) 

which has been proposed as a good candidate to describe the evolution of the wave function of 
Josephson junctions 1 79 ]. The corresponding Weyl symbol is 

h({,4>) = — (— +Acos<f) . (69) 



2n\2 

Using the relations ( |60) , we get the evolution equation for the Wigner function 

d 



-2^sin^>sinh((5 f /2) 
2 d(p 



W(€,</>). (70) 



We can directly check that exp(6 M /2)W(f, 0) = W(€ + 1/2, (p), so |70| becomes 

| W, <f>) = ~ l~ W(f, <f>)-A sin ftWtf + 1/2,0) - W(€ - 1/2, 0)] ] . ( 7 1 ) 



We are interested in the semiclassical evolution, i.e., for states whose angular momentum 
components are sufficiently concentrated around a certain value (q with » 1. Thus, we may 
view W{£ + 1/2, </>) - W(( - 1/2, <p) as a derivative respect to t, so we finally get 



i^-^d^-^i)^' (72) 

whose formal solution is W(i( f), <f>(t)) in terms of the "classical" trajectories 

^ = -2i sin 0, (f> = -£. (73) 
This constitutes a nice solution for an involved problem. 

4. 1 . Tomography 

To complete our theory we propose a reconstruction scheme for observables in this cylindrical 
phase space. The equivalent version of pT) reads now, particularized for the density operator, 



=2 [V,< 

pel? 



) £>(€, <f>) d<f> , (74) 



where q({, </>) = Tt[q ffi(€, </>)]/(2n). In terms of q{(, <f>), the Wigner function is 

W,0 = ;r£ f expW€>-^')]^',0')#'. (75) 

A reconstruction of this Wigner function is thus tantamount to finding the coefficients q((, <p). 
For £ = 0, these coefficients 

e(0, 0) = 2 «p(^) < W ( ?6 > 

71 eeZ 

are simply the Fourier transform of the angular-momentum spectrum. 

For t ± 0, we assume that we are able to experimentally determine L 2 transformations on the 
input state followed by angular projections, that is, 

p(<P, = <0| e X p(/^L 2 /2)£ exp(-#£ 2 |0> . (77) 

One can check that 

e(A 0) = Y n \ e*&rW)ptf* <PI ( ) , (78) 
so the measurement of p((p, Q allows for the determination g{{, <p) and hence the full reconstruc- 



tion of the Wigner function via equation ( 75 1 



As a rather simple yet illustrative example, let us note that for the vortex state |€o), L is 
diagonal, so the tomograms p{$ ,<$>' j€) are independent of <p an d t an d all of them equal to 
1/(2tt). Performing the integration we obtain precisely the Wigner function in equation ( |6T| >. 

The feasibility of the proposed scheme relies on two crucial points. First, the implementation 
of the L 2 transformation, which corresponds to a free rotor. This has been used to describe the 
evolution in a variety of situations. Second, we need to assess the measurement of the angular 
spectrum, which can be done only approximately. Though the implementation of this scheme 
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may differ depending on the system under considerations, our formulation provides a common 
theoretical framework on the Hilbert space generated by the action of angle and angular momen- 
tum. A experimental demonstration of the method in terms of optical beams has been recently 
established ITHOI ■ 

For continuous variables, it is possible to make a direct sampling of the Wigner function by 
a scheme that determines the parity [81 82 1. This is an unexplored territory for the cylinder that 
is worth investigating. 

5. Concluding remarks 

In summary, we have shown how to extend in a consistent way all the techniques developed 
for a continuous-variable phase space to the case of angle and angular momentum. While we 
have not left aside the mathematical details, our main emphasis has been on presenting a simple 
and useful toolkit that any practitioner in the field should master. In our view, far from being 
an academic curiosity, the ideas expressed here have a wide range of potential applications in 
numerous hot topics. 

This paper has greatly benefited from the criticism and advice of Prof. B.-G. Englert. The 
work was supported by the Spanish Research Directorate (Grants FIS2005-06714 and FIS2008- 
04356), the UCM-BCSH program (Grant GR-920992), the Mexican CONACyT (Grant 106525), 
the Czech Ministry of Education (Project MSM6 1989592 13), and the Czech Ministry of Industry 
and Trade (Project FR-TI1/364). 

References 

[1] H. Weyl, Gruppentheorie und Quantemechanik, Hirzel-Verlag, Leipzig, 1928. 

[2] E. P. Wigner, Phys. Rev. 40(1932)749-759. 

[3] J. E. Moyal, Proc. Cambridge Phil. Soc. 45 (1949) 99-124. 

[4] R. L. Stratonovich, JETP 31 (1956) 1012—1020. 

[5] G. S. Agarwal, E. Wolf, Phys. Rev. D 2 (1970) 2161-2186. 

[6] F. A. Berezin, Commun. Math. Phys. 40 (1975) 153-174. 

[7] G. S. Agarwal, Phys. Rev. A 24 (1981) 2889-2896. 

[8] J. Bertrand, P. Bertrand, Found. Phys. 17 (1987) 397^105. 

[9] J. C. Varilly, J. M. Gracia-Bondfa, Ann. Phys. 190 (1989) 107-148. 

[10] C. Brif, A. Mann, J. Phys. A 31 (1998) L9-L17. 

[11] M. G. Benedict, A. Czirjak, Phys. Rev. A 60 (1999) 4034-^1044. 

[12] V. I. Tatarskii, Sov. Phys. Usp. 26 (1983) 31 1-327. 

[13] N. L. Balazs, B. K. Jennings, Phys. Rep. 104 (1984) 347-391. 

[14] M. Hillery, R. F. O. Connell, M. O. Scully, E. P. Wigner, Phys. Rep. 106 (1984) 121-167. 

[15] H.-W. Lee, Phys. Rep. 259 (1995) 147-211. 

[16] F. E. Schroek, Quantum Mechanics on Phase Space, Kluwer, Dordrecht, 1996. 

[17] W. P. Schleich, Quantum Optics in Phase Space, Wiley-VCH, Berlin, 2001. 

[18] C. K. Zachos, D. B. Fairlie, T. L. Curtright (Eds.), Quantum Mechanics in Phase Space, World Scientific, Singapore, 
2005. 

[19] W. K. Wootters, Ann. Phys. 176 (1987) 1-21. 

[20] D. Galetti, A. F. R. de Toledo-Piza, Physica A 149 (1988) 267-282. 

[21] D. Galetti, A. F. R. de Toledo-Piza, Physica A 186 (1992) 513-523. 

[22] W. K. Wootters, IBM J. Res. Dev. 48 (2004) 99-1 10. 

[23] K. S. Gibbons, M. J. Hoffman, W. K. Wootters, Phys. Rev. A 70 (2004) 062101 . 

[24] A. Vourdas, Rep. Prog. Phys. 67 (2004) 267-320. 

[25] A. B. Klimov, C. Munoz, J. L. Romero, J. Phys. A 39 (2006) 14471-14497. 

[26] A. Vourdas, J. Phys. A 40 (2007) R285-R33 1 . 

[27] G. Bjork, A. B. Klimov, L. L. Sanchez-Soto, Prog. Opt. 51 (2008) 469-516. 



14 



[28] A. B. Klimov, J. L. Romero, G. Bjork, L. L. Sanchez-Soto, Ann. Phys. 324 (2009) 53-72. 

[29] T. Durt, B.-G. Englert, I. Bengtsson, K. Zyczkowski, Int. J. Quantum Inf. 8 (2010) 535-640. 

[30] M. G. G. Laidlaw, C. M. DeWitt, Phys. Rev. D 3 (1971) 1375-1378. 

[31] L. S. Schulman, J. Math. Phys. 12 (1971) 304-308. 

[32] J. S. Dowker, J. Phys. A 5 (1972) 936-943. 

[33] F. Acerbi, G. Morchio, F. Strocchi, Lett. Math. Phys. 27 (1993) 1-11. 

[34] L. Allen, S. M. Barnett, M. J. Padgett, Optical Angular Momentum, Institute of Physics Publishing, Bristol, 2003. 

[35] D. L. Andrews (Ed.), Structured Light and Its Applications, Academic Press, Boston, 2008. 

[36] G. Molina-Terriza, J. P. Torres, L. Torner, Nat. Phys. 3 (2007) 305-310. 

[37] S. Franke- Arnold, L. Allen, M. Padgett, Laser Photon. Rev. 2 (2008) 299-313. 

[38] A. Vaziri, G. Weihs, A. Zeilinger, J. Opt. B 4 (2002) S47-S51. 

[39] N. K. Langford, R. B. Dalton, M. D. Harvey, J. L. O'Brien, G. J. Pryde, A. Gilchrist, S. D. Bartlett, A. G. White, 

Phys. Rev. Lett. 93 (2004) 053601. 

[40] P. Carruthers, M. M. Nieto, Rev. Mod. Phys. 40 (1968) 41 1-440. 

[41] R. Lynch, Phys. Rep. 256 (1995) 367^36. 

[42] V. Pefinova, A. Luks, J. Pefina, Phase in Optics, World Scientific, Singapore, 1998. 

[43] A. Luis, L. L. Sanchez-Soto, Prog. Opt. 44 (2000) 42 1-48 1 . 

[44] N. Mukunda, Am. J. Phys. 47 (1979) 182-187. 

[45] N. Mukunda, G. Marmo, A. Zampini, S. Chaturvedi, R. Simon, J. Math. Phys. 46 (2005) 012106. 

[46] J. P. Bizarro, Phys. Rev. A 49 (1994) 3255-3276. 

[47] A. Vourdas, J. Phys. A 29 (1996) 4275^1288. 

[48] L. M. Nieto, N. M. Atakishiyev, S. M. Chumakov, K. B. Wolf, J. Phys. A 31 (1998) 3875—3895. 

[49] M. Ruzzi, D. Galetti, J. Phys. A 35 (2002) 4633^640. 

[50] S. Zhang, A. Vourdas, J. Math. Phys. 44 (2003) 5084-5094. 

[51] K. Kakazu, E. Sakai, Prog. Theor. Phys. 115 (2006) 1027-1045. 

[52] K. Kowalski, J. Rembielihski, L. C. Papaloucas, J. Phys. A 29 (1996) 4149^167. 

[53] J. A. Gonzalez, M. A. del Olmo, J. Phys. A 31 (1998) 8841-8857. 

[54] Y. Ohnuki, S. Kitakado, J. Math. Phys. 34 (1993) 2827-2851. 

[55] B. C. Hall, J. J. Mitchell, J. Math. Phys. 43 (2002) 121 1-1236. 

[56] M. Ruzzi, M. A. Marchiolli, E. C. da Silva, D. Galetti, J. Phys. A 39 (2006) 9881-9890. 

[57] I. Rigas, L. L. Sanchez-Soto, A. B. Klimov, J. Rehacek, Z. Hradil, Opt. Spectrosc. 108 (2010) 206-212. 

[58] H. A. Kastrup, Phys. Rev. A 73 (2006) 052104. 

[59] I. Rigas, L. L. Sanchez-Soto, A. B. Klimov, J. Rehacek, Z. Hradil, Phys. Rev. A 78 (2008) 060101. 

[60] E. Binz, S. Pods, The Geometry of Heisenberg Groups, American Mathematical Society, Providence, 2008. 

[61] A. Galindo, P. Pascual, Quantum Mechanics, Springer, Berlin, 1991. 

[62] A. Kenfack, K. Zyczkowski, J. Opt. B 6 (2004) 396-404. 

[63] R. L. Hudson, Rep. Math. Phys. 6 (1974) 249-252. 

[64] A. J. E. M. lanssen, SIAM J. Math. Anal. 15 (1984) 170-176. 

[65] E. H. Lieb, J. Math. Phys. 31 (1990) 594-599. 

[66] A. Perelomov, Generalized Coherent States and their Applications, Springer, Berlin, 1986. 

[67] M. L. Mehta, J. Math. Phys. 28 (1987) 781-785. 

[68] M. Ruzzi, J. Math. Phys. 47 (2006) 063507. 

[69] M. G. A. Paris, J. Rehacek (Eds.), Quantum State Estimation, volume 649 of Led. Not. Phys. , Springer, Heidelberg, 
2004. 

[70] A. I. Lvovsky, M. G. Raymer, Rev. Mod. Phys. 81 (2009) 299-332. 

[71] G. G. Emch, Algebraic Methods in Statistical Mechanics and Quantum Field Theory, Wiley, New York, 1972. 

[72] J. B. M. Uffink, Measures of Uncertainty and the Uncertainty Principle, Ph.D. thesis, University of Utrecht, 1990. 

[73] R. Bracewell, The Fourier Transform and Its Applications, McGraw-Hill, 3rd ed edition, 1999. 

[74] Z. Hradil, J. Rehacek, Z. Bouchal, R. Celechovsky, L. L. Sanchez-Soto, Phys. Rev. Lett. 97 (2006) 243601. 

[75] J. Rehacek, Z. Bouchal, R. Celechovsky, Z. Hradil, L. L. Sanchez-Soto, Phys. Rev. A 77 (2008) 0321 10. 

[76] J. F. Plebahski, M. Prazanowski, J. Tosiek, F. K. Turrubiates, Acta Phys. Pol. B 31 (2000) 561-587. 

[77] D. Mumford, Tata Lectures on Theta I, Birkhauser, Boston, 1983. 

[78] I. Rigas, L. L. Sanchez-Soto, A. B. Klimov, J. Rehacek, Z. Hradil, Phys. Rev. A (2009). 

[79] J. R. Anglin, P. Drummond, A. Smerzi, Phys. Rev. A 64 (2001) 063605. 

[80] J. Rehacek, Z. Hradil, Z. Bouchal, A. B. Klimov, I. Rigas, L. L. Sanchez-Soto, Opt. Lett. 35 (2010) 2064-2066. 

[81] B.-G. Englert, N. Sterpi, H. Walther, Opt. Commun. 100 (1993) 526-535. 

[82] P. Bertet, A. Auffeves, P. Maioli, S. Osnaghi, T. Meunier, M. Brune, J. M. Raimond, S. Haroche, Physical Review 
Letters 89 (2002) 200402. 



15 



